Searching for Non-Abelian Phases in the Bose-Einstein Condensate of Dysprosium 
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The recently realized Bose-Einstein condensate of Dysprosium will become a spin-8 spinor con- 
densate at ultralow magnetic fields. In such a high spin condensate, many phases with different 
symmetries can exist. Among them the most interesting ones are those with non-abelian point 
group symmetry. In this letter we discuss the variety of symmetry phases in a spin-8 condensate 
resulting from numerical solutions of the Hamiltonian. We show that these symmetries can be 
determined uniquely from the measurements of density population on each spin component in an 
ultralow magnetic field, together with the measurements of the collective modes in the zero-field 
limit. This method can also be applied to Bose-Einstein condensate of other magnetic atoms, such 
as Cr and Er. 



Recently Benjamin Lev's group has successfully re- 
alized Bose-Einstein condensate of Dysprosium (Dy)[T]. 
Dy is a complex lanthanide atom with totally 10 /-shell 
electrons. It has electron spin 5 = 2, orbital angular mo- 
mentum L = 6, and nuclear spin 1 = 0. The total spin 
F = L + S becomes F = 8 due to the Hund's rule. The 
spin of Dy can be easily polarized by an external mag- 
netic field, which will then give rise to strong magnetic 
dipolar interactions. However, it is technically possible to 
reduce the magnetic fields to very low values so that the 
system is essentially depolarized, as recently achieved in 
the Cr condensate [2 . In this limit, the spin Hamiltonian 
has rotational and gauge symmetry, G = SO (3) x U(l). 
Bose condensation breaks this symmetry. However, the 
resulting state can have the symmetry of a subgroup of 
G, which can be non-abelian for large spin bosons [3]. 
Condensates with non-abelian symmetry group are par- 
ticularly interesting, as they have non-abelian topological 
defects with novel physical properties [4 . 

At zero- field limit, the actual spin ground state of Dy 
is determined by nine scattering lengths that are not yet 
known. It will be overcomplicated to construct a phase 
diagram with all these unknown parameters. In this let- 
ter, we focus on the question of non-abelian phases of 
Dy condensate, and the simplest scheme to detect them. 
While the non-abelian nature of the order parameter is 
fully manifested in the resulting structure after its line 
defects cross each other [4 , such experiments are too in- 
volved in a cold atom setup at present. In this work, we 
propose a much simpler method to determine the sym- 
metry of a spinor condensate. It combines the measure- 
ments of spin populations and the Bogoliubov spectrum. 
The former can be achieved with Stern- Gerlach technique 
while the later can be measured by Bragg spectroscopy. 
Both techniques are widely used in cold atom experi- 
ments today. 

Model: Let i/j m be the bosonic operator for each spin 
component m, (m = —F, . . . , F), the Hamiltonian for Dy 



(F = 8) is = tfo + tfint [3 El, 
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where gj = 27th 2 a j /M ', aj is the s-wave scattering length 
between a pair of bosons with total spin j, described by 
the local pair operator A/ m (r) = ^2 m (j, m|F, mi, F, m — 
mi)^ TOl (r)^ TO _ TOl (r). Bose statistics, however, implies 
j = 0,2,. .16 only. B is a tiny magnetic field along z 
direction, and the magnetic moment \i = 5/ie/4 for Dy 

The energy of a condensate (VVn,k=o) = ^/no^Pm is 

2F F 

£= ^2 gjnl^2\A jrn \ 2 -n Bfi ^ m|(/? m | 2 (3) 
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where Aj m (j, m|F, mi, F, m - mi)<p 

while no represents the total density. Using the imagi- 
nary time evolution method, we have conducted exten- 
sive searches of possible ground states in zero and small 
fields for vast ranges of gj (j = 0, 2, . . . , 2F). We have fur- 
ther studied the Bogoliubov spectra of the ground states 
to ensure their stability, which also turns out to be very 
useful for determining the symmetry of the state. 

Majorana Representation: A very useful way to de- 
scribe the condensate wavefunction {tp m } is to use the 
Majorana representation of spin states [8]. Recently, 
many authors have applied this method to study Bose 
condensate with high spins [3j EJ [10]. This repre- 
sentation is most conveniently described in terms of 
Schwinger bosons [10 j. This method presents the spin 
operators in terms of bosons a and 6, such that F x = 
(a)b + 6 t a)/2, F y = i(a)b - 6 t a)/2 and F z = (a) a - 
b^b)/2. A general normalized spin state can be written as 
W) = T, F m=- F ^m\F,m), where \F,m) = [(F + m)\(F- 
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FIG. 1. Major ana representations of spin-8 states with dif- 
ferent symmetries. Their descriptions are given in the dis- 
cussion of Symmetry Allowed States. The point with po- 
lar angle (9i,<fii) on the unit sphere represents the spinor 
( Ui ,Vi) = (cos^/2,e^sin^/2). 

m)\]- 1 / 2 a^ F+m ^ F - m \0). Since the sum in \<p) is a ho- 
mogenous polynomial of a) and £r of degree 2F, it can 
be factorized as 

F 1 2F 

\<p) = Vrn\F,m) = +^ f )l°)' ( 4 ) 

771= ~F 2=1 

where Af is a normalization factor to ensure all \ui\ 2 + 
|^| 2 = 1. Absorbing the overall phase factor of the 
product into A/", we can write Ui — cos(# z /2) and v\ = 
e 1 ^ 1 sin(^/2), so that the spinor (ui,Vi) can be repre- 
sented as a point on a unit sphere denoted by a unit 
vector with polar angle (0^,(^). It is straightfor- 
ward to show that (i) under a spin rotation, n 2 - ro- 
tates as a Cartestian vector, and (ii) under time rever- 
sal, (iii,Vi) — > (v*j— u*), which implies n 2 - = — n 2 . The 
2F + 1 complex numbers {(/? m } (with totally AF + 1 real 
variables because of normalization) is now replaced by 
2F unit vectors n 2 - plus an overall phase factor. 

Symmetry Allowed States: If U is a rotational sym- 
metry operation, it means the state returns to itself up 
to an overall phase after being operated on by J7, i.e. 
-> Y,m> U mm>(Pm' = e za (p m . Since the symmetry 



group G of H is 50(3) x £7(1), the symmetry of the 
ground state will be described by one of the normal point 
groups E, O n , D ni T, O or F[4]. E 1 is the identity, corre- 
sponding to states with no particular symmetry. O n de- 
notes an n-fold rotational symmetry along certain axis. 
D n contains one n-fold C n axis plus n two-fold O2 axes 
perpendicular to O n . T, O, and Y denote tetrahedral, 
octahedral, and isocahedral symmetry respectively. The 
following properties are useful for our later discussions. 

(I) If the vectors {n z } comes in (n, — n) pairs, the state 
is invariant under time reversal. 

(II) A ferromagnetic state (by which we mean (F) ^ 0) 
breaks time reversal symmetry and can have only one 
O n axis. States with more than one C n axis such as 
those with D n , T, O, and Y symmetry can not be fer- 
romagnetic, even though they may break time reversal 
symmetry. 

Now let us analyze all the symmetry states: 

(1) E: These states have no special symmetry. They 
are not found as ground states in our energy minimiza- 
tion. 

(2) O n : Dy has F = 8, hence 16 points on the 
unit sphere. So we either have n ^ 16 or n = 00. The 
O n groups can be further classified into C nvi C n h and 
52n[H]- We have found several examples of C nv symme- 
try. For instance, a C$ v state is shown in Fig. [TJa), in 
which all 16 points form three pentagons and one point 
sitting at the south pole. However, examples with only 

symmetry are not found [12 . We also get 
states with symmetry, with q points (q < 8) collaps- 
ing into one and other 16 — q points into the antipodal 
point. The case of symmetry with q = 1 is shown 
in Fig. [TJb). All the states with C nv and symmetry 
break time reversal symmetry, and have ferromagnetic 
order. 

(3) D n : In this case, either n = 00 or n is a finite 
integer such that n ^ 16 for even n, and n ^ 7 for odd 
n [13] . D n symmetry group can be further classified into 
D n d and D n h[ll r In Fig. [TJc) and (d), we show two 
examples of D n &. D$d consists of two octagons, one in the 
north hemisphere and the other in the south hemisphere. 
D^d consists of four squares, two in the north hemisphere 
and the other two in the south hemisphere. In Fig. [lje), 
we show an example of F) nh . For -Du^, 14 points are 
evenly distributed on the equator and the other two are 
placed on the north and south poles. While some of 
these states break time reversal symmetry, such as D n d, 
none of the states with D n symmetry are ferromagnetic. 
We, however, do not find the "polar" state with 8 points 
collapse into the north and the other 8 points to the south 
pole. We suspect it either does not exist or only occupies 
a very limited region in the phase diagram. 

(4) T: We have found one phase with T-symmetry and 
another with Td symmetry[ll j. They occupy a large por- 
tion of the phase diagram. An example of the state with 
Td symmetry is shown in Fig. [ljf ) . Four points are dis- 
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tributed at the vertices of a tetrahedron. Around each 
tetrahedral vertex, another three points are distributed 
at equal distance d from it forming a small regular trian- 
gle. The three vertices of this triangle lie in the three 
mirror planes passing through this tetrahedral vertex. 
The distance d varies with interaction and can shrink 
to zero, as shown in Fig. [I] (g). We have also found a 
phase with T symmetry where the vertices of the trian- 
gles move away from the mirror plane, as shown in Fig. 
[TJh). This transition from to T is second order. Both 
phases T and Td break time reversal symmetry. 

(5) O: This state is described by a cube with two points 
at each Conner, as shown in Fig. [IJi). This state does 
not break time reversal symmetry, and is therefore non- 
magnetic. 

(6) Y: This is impossible for Dy, as one can not dis- 
tribute 16 points on a sphere with icosahedral symmetry. 
On the other hand, such a phase is possible for Erbium 
(Er) which is F = 6 and is represented by 12 points on 
the sphere [T4]. 

All the states with D n (n finite), T or O symmetry have 
non-abelian defects [4]. Hereafter we will discuss how to 
detect them experimentally. 

Detecting the Symmetry: We propose two measure- 
ments (A) and (B) for probe the symmetry of the system. 
Measurement (A) is to pin a C n axis by an ultralow mag- 
netic field along z. We find numerically that such pinning 
always happens [7 . The spin wavefunction along z must 
then be of the form 

<p =(.-., 0, a, Or^O, ft (V^A 7, 0, • • • ) T (5) 

n — 1 n — 1 

for any finite n, i.e. any non- vanishing component must 
be separated from the next non- vanishing one by n — 1 ze- 
ros. If n = oo, there will be only one non-zero component. 
The positions and the magnitudes of the non-vanishing 
components are also different for different states. This 
follows from the fact that a spin rotation along z axes by 
2ir/n will change ip to e %e (p, hence e~* 2m7r / n (/? m = e 2(9 (/? m , 
which can only be satisfied if ip takes the form of Eq. (J5|. 

Measurement (B) concerns the low-energy collective 
modes in zero-field limit. As explained later, for nonmag- 
netic states with discrete point group symmetry, there 
will be totally four gapless Goldstone modes with linear 
dispersion. One corresponds to the broken global U(l) 
symmetry, while the other three correspond to broken 
spin rotation symmetry. If the remaining symmetry con- 
tains only one C n> 2 axis, two spin rotational modes will 
be degenerate (i.e. only three different dispersions). If 
there are more than one C n> 2 axis, all three rotational 
modes will be degenerate (i.e. two different dispersions). 
For ferromagnetic ground states, there will be a quadratic 
Goldstone mode, as already known from previous study 
of magnetism [5 . While these results are found from ex- 
plicit calculations of the Bogoliubov modes, they in fact 
result from the following general considerations. 




FIG. 2. Low-energy spectra of Bose condensates with different 
symmetries, k is in unit of ko = y/Tmjl/h and uj is in unit of 
0.1/i (/jl: chemical potential). The green (red) line stands for 
two (three) degenerate modes. See also Tab. u\ 



In the low-energy limit k — >• 0, the three spin Gold- 
stone modes correspond to simultaneous rotation of all 
vectors along axes x, y and z. If there is one C n> 2 
axis, say, along z, such a rotation transforms the fluctua- 
tions along x and y to those along x' and y' \ respectively. 
Due to the C n symmetry, the fluctuations along x (y) 
will behave the same as those along x' (y f ). On the other 
hand, the fluctuations along x' or y' are linear combina- 
tions of those along x and y. It then implies these two 
modes along x and y must be degenerate. Furthermore, 
if there is another C n> 2 axis, say, along x direction, it 
will imply that fluctuations along y and z also have the 
same dispersion. Then, all three rotational modes are de- 
generate. These results can be applied to different point 
groups as summarized in Tab. [I| 

Experimental Scheme for Differentiating Different 
States: Applying the Methods (A) and (B), we can make 
the following distinctions. 

(i) Distinguishing E and C n from others: As seen in 
Tab. |TJ These are the only two ferromagnetic states, and 



TABLE I. Illustration of collective modes of different symme- 
try states. "L" and "Q" stand for modes linear and quadratic 
in k respectively. "2(3) dg." is short for two (three) degener- 
ate modes. 



Symmetry 


Phase Mode 


Spin Modes 


Distinct Dispersions 


E 


1L 


1L + 1Q 


2L + 1Q 


C n 


1L 


1L + 1Q 


2L + 1Q 


D n 


1L 


3L (2 dg.) 


3L 


T,0 


1L 


3L (3 dg.) 


2L 


Doo 


1L 


2L (2 dg.) 


2L 
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One mode 
quadratic 



All modes 
linear 



3 different 
velocities 



All \ (p m \ 2 non-zero 




n-1 in every n consecutive \(p m \ 2 vanish 



D„ 



2 different 
velocities 



Only one|<p m | 2 ^ or 1 in 2 
consecutive | (p m | 2 vanishes 



2 in 3 consecutive | (p m | 2 vanish 



3 in 4 consecutive | (p m | 2 vanish 



T 




Measure Bogoliubov mode: 



Measure spin population 



FIG. 3. Schemes for determining the symmetry of different 
phases. Blue and red arrows denote the measurement of col- 
lective modes and spin populations \(f m \ 2 respectively. 



thus have quadratic Goldstone modes. The case of C$ v 
is show in Fig. [2|a). Thus, Measurement (B) can distin- 
guish them from other states. Moreover, E and C n can 
be distinguished and determined by Measurement (A), 
as C n has the characteristic order parameter as shown in 
Eq.(j5|. (Note that our calculations show that an external 
magnetic field always pin the C n axis of C n states). 

(ii) Distinguishing D n from D^, T and O: As seen 
in Tab. [I| D n differs from T and O from the number 
of mode velocities. (See also Fig. J2jb)-(d)). In mea- 
surement (A), the untralow magnetic field pins along z 
either C n or C2 axis of the D n symmetry, depending on 
the states [7i . If C2 axis is pinned, then we can add a 
small magnetic field gradient B = B[z + Go(xx — zz)\. 
This is effectively equivalent to adding an energy e(i^), 
with e = h 2 Gl/(2M) [15], which will pin the C n axis 
along either x or y direction. One can then determine n 
by measuring the spin population along x or y [7] . 

(hi) Distinguishing .Doo, T and O: Our calculations 
show an ultralow magnetic field always pin the C3 axis of 
T symmetry, while pins the C4 axis of the O state. Even 
though our calculations did not find the "polar" state 
with Dqq symmetry, we note that if present, it can only 
be pinned along either the or the C2 axis by magnetic 
field, which is distinct from T and O. So by measurement 
(A) we can easily distinguish these three. Another dif- 
ference between T and O is that T phase breaks time 
reversal symmetry, while O phase does not. Hence for O 
phase, we have |(/? m | 2 = |(/?_ m | 2 in zero field limit, which 
is a property not shared by the T phase [7] . 

In Fig. [3] we summarize our proposal with a flow di- 
agram. Following these steps, all possible (normal point 
group) symmetries can be determined exactly. 

Final Comments: Although our discussions are for Dy, 
our method of detection is applicable for all spinor con- 
densates such as Cr [2], and Er [14]. In our discussion, 
we have not considered dipole energy. The competition 
between dipole energy and spin-dependent s-wave inter- 
action will certainly change the energetics. However, 
dipolar energy is known to depend strongly on the ge- 



ometry of the trap, and is known to have less effects in 
spherical potentials. Moreover, since all the non-abelian 
states are nonmagnetic, dipolar interaction may not have 
strong effects on them. In settings where dipolar inter- 
action becomes dominant, the system tends to develop 
non-uniform spin textures as shown in the recent exper- 
iments on 87 Kb [T6] . Even in such cases, the underlying 
zero field states such as those discussed here still play 
an important role in determining the global equilibrium 
spin textures [T7j. 
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SUPPLEMENTARY MATERIAL 
Magnetic Moment of Dysprosium 

Dysprosium (Dy) atom has 10 / shell electrons, which give the atom electron spin S = 2 and electron orbital 
angular momentum L = 6 according to the first and second Hund's rule. Due to spin orbit coupling, the total spin 
F = L + S of Dy atom becomes F = 8 (the third Hund's rule). The magnetic moment operator of the atom is then 
/x = /i&(L + 2S), where /i& = 5.79 x 10 _5 eV • T _1 is Bohr magneton. Typically the spin orbit coupling energy of an 
atom is around 10 _2 eV, while Zeeman energy and interaction energy per atom in experiments are usually lower than 
10 _8 eV. Thus Dy atoms will nearly perfectly stay in the coupled F — 8 representation, and we can use the projection 
of /i into the F — 8 subspace as the magnetic moment. 

Both F and {i are SU(2) irreducible tensor operators of order 1 (vector operators). Wigner-Eckart theorem says an 
order k SU (2) irreducible tensor operator T£ behaves as 

{F',rn^\F,rn) = ^^^{F'\\V^F) (6) 

where {F'm'kq\Fm) is the Clebsch-Gordan coefficient, and {F'\ \ V k \ \F) is a factor independent of 777/ , m, q. For F' = F, 
Eq. (|6| then tells us all the tensor operators of spin k are proportional when projected into the spin F subspace. As 
a special case of k = 1, the projected magnetic moment [i is proportional to spin F, with a factor \i to be determined 
below. 

By a dot product with F on both sides of fi = /iF, we arrive at the following equation 

Mb (L + S) • (L + 2S) = [^F 2 - i(L 2 - S 2 )] = /jF 2 (7) 
By noticing L 2 = L(L + 1), S 2 = S(S + 1) and F 2 = F(F + 1), we can derive out 

^-^2 2F(F + 1) ) (8) 

For Dy atom with S = 2, L = 6 and F = 8, we have \i = 5/i^/4. 



Interaction Parameters and Ground State 

By using the relationship Y^j=o 2 ••• Ei=-j l^j'm| 2 = (^V) 2 ? ^ ne mean- field interaction energy can be rewritten as 

2F 

£int = ^O^o(^V) 2 + Yl ^°) n Yl I A'm| 2 (9) 

j=0,2,"- m 

where A jm E mi {j,rn\F,m u F,m - mi)<p m 1 ( Pm—m 1 - Only the second part of Eq. is spin dependent and useful 
for the determination of the ground state, which consists of 8 parameters gj — go Tor F = 8. As examples, the 
interaction parameters which lead to ground states shown in Fig. 1 of our paper are summarized in Table |TT| (without 
external magnetic field). These parameters {gj — go} are chosen to be located at the interior of the corresponding 
phases' space, and therefore some variation of these parameters will not change the symmetry of the states. One 
should also note that a simultaneous multiplication on {gj — go} by a factor A (i.e. {gj — go} — » {A(#j — go)}) do not 
change the ground state. 
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Spin Population in a Magnetic Field 

When an ultralow magnetic field B = Bz is turned on, the symmetry of Hamiltonian is reduced to SO (2) x U(l), 
and the ground state will generally have one of its rotational axis fixed in z direction. For C n symmetry which 
are ferromagnetic phases, the spin (F) points along C n axis. Thus the first order perturbation energy — /iB • (F) 
is non-zero, and the C n axis will be fixed along z direction. For non-magnetic phases with D ni T or O symmetry, 



it is the second order perturbation that determines which axis is fixed toward z. In Table III the calculated spin 
populations (under z basis) of several phases in a ultralow magnetic field are shown, from which it is easy to figure out 
the fixed axis according to Measurement (A). In the calculations the Zeeman energy added is around 1% ~ 10% of 
the spin dependent interaction energy, and we note this Zeeman energy only changes the magnitude of non-zero spin 
components, while keeps vanishing components (in the zero field limit) zero. For different states with D n symmetry, it 
is either the C n or the C2 axis that is fixed along z direction. Numerical computations show for D n d (D nv ) symmetry, 
C2 (C n ) axis is more likely fixed. For T and O symmetry, C3 axis and C4 axis are fixed along z, respectively. 



Effect of Magnetic Field Gradient for D nd Symmetry 

For D n phase whose C2 axis is fixed along z (such phases are likely to have a D n d symmetry), measuring the spin 
populations under z basis does not give us the number n. To figure out n, we need to add an ultralow magnetic field 
with small gradient B = B[z + Go(xx — zz)\. By a local unitary transformation (p — >• e~ lG ° FyX p, all local z axes are 
aligned along B, and a kinetic energy cost [15] 

S k = -n h 2 ^e iGoF y x V 2 e- iGoF y x cp/(2M) = n e(F 2 ) (10) 

arises, where e = h 2 Gl/(2M). This kinetic energy must be far smaller than Zeeman energy, for otherwise the Bose 
Einstein Condensate will prefer the uniform configuration that costs the Zeeman energy instead of the kinetic energy. 
This kinetic energy term will fix C n axis of the D n symmetry along either x or y. First we know for D n phases with 
one of the C2 axes fixed along z by external field B, the C n axis then lies in x — y plane. By symmetry consideration, 
C n must be one principal axis of the real symmetric tensor = ((FiFj+FjFi)} /2 (i, j = x, z). Thus (F 2 ) = y-Q-y 
will reach the minimum when C n axis points along either x or y. A Measurement of the spin populations in either x or 
y basis will then give us the number n. For the two phases and shown in Fig. 1(c) and (d), our calculations 
show that their C$ and C4 axes are fixed along x and y, respectively. 

In experiments it is hard to change the basis for spin population measurement, so an easier method is to perform a 
7r/2 rotation of atoms about y or x axis, followed by a measurement of spin population in z basis. A tt/2 rotation about 
y axis can be realized by adding an additional alternating magnetic field B a = B a y cos Qt where Q is the Zeeman 
splitting frequency fiB/H. The frequency Q will equivalently recover the degeneracy of split Zeeman levels, and the 
atoms will rotate an angle fiB a T/h about y axis within time T. tt/2 rotation is realized by controlling T = ttH/ (2fiB a ). 

I 



TABLE II. Interaction parameters g^j — go which lead to ground states shown in Fig. 1. 



Phase 


Isotropy group 


#2 - go 


#4 - go 


g& - go 


#8 - go 


010 - go 


gi2 - go 


0i4 - go 


016 - go 


Fig. 1(a) 




1 








-1 


-2 





-2 


1 


Fig. 1(b) 


Coo 


1 








1 


4 





-5 


1 


Fig. 1(c) 


D 8d 


-1 




















1 


Fig. 1(d) 


D 4d 


1 








1 


4 








1 


Fig. 1(e) 




2 


2.5 





1.5 


5 


7 


2.5 


4 


Fig. 1(f) 


T d 


0.1 








-1 


-1 








1 


Fig. 1(g) 


T d 


0.1 








-1.27 


-1 








1 


Fig. 1(h) 


T 


0.1 








-8 


-1 








1 


Fig. l(i) 


o h 


1 


1 


1 


1 


-5 


-5 


1 


1 
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TABLE III. Calculated population p m — tptnPm measured along -\-z direction weak magnetic field. In the calculations we are 
adding a Zeeman energy around 1 ~ 10 percent s of the spin dependent interaction energy. 



Phase 


Fig. 1(a) 


Fig. 1(b) 


Fig. 1(c) 


Fig. 1(d) 


Fig. 1(e) 


Fig. 1(f) 


Fig. l(i) 


Isotropy group 




Cooy 


D 8d 




Di4h 


T d 


O 


Axis fixed 


c 5 


Coo 


c 2 


c 2 


Cl4 


C 3 


c 4 


P8 








0.1191 


0.1963 








0.2210 


P7 


0.4922 


1.0000 








0.5590 


0.0003 





P6 








0.0649 


0.4341 











P5 























P4 








0.0780 


0.0170 





0.4953 


0.0717 


P3 























P2 


0.4629 





0.1559 


0.2254 











Pi 

















0.1736 





Po 








0.2110 


0.0583 








0.5179 


P-l 























P-2 








0.1451 


0.0444 





0.0604 





P-3 


0.0424 




















p-4 








0.0693 


0.0005 








0.0646 


P-5 

















0.0426 





p-6 








0.0567 


0.0192 











P-7 














0.4410 








P-8 


0.0025 





0.1000 


0.0049 





0.2279 


0.1247 



